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ABSTRACT 

It lias been posited that lunar eclipse observations may help predict the in-transit 
signature of Earth-like extrasolar planets. However, a comparative analysis of the two 
phenomena addressing in detail the transport of stellar light through the planet's at- 
mosphere has not yet been presented. Here, we proceed with the investigation of both 
phenomena by making use of a common formulation. Our starting point is a set of pre- 
viously unpublished near-infrared spectra collected at various phases during the August 
2008 lunar eclipse. We then take the formulation to the limit of an infinitely distant 
observer in order to investigate the in-transit signature of the Earth-Sun system as be- 
ing observed from outside our Solar System. The refraction-bending of sunlight rays 
that pass through the Earth's atmosphere is a critical factor in the illumination of the 
eclipsed Moon. Likewise, refraction will have an impact on the in-transit transmission 
spectrum for specific planet-star systems depending on the refractive properties of the 
planet's atmosphere, the stellar size and the planet's orbital distance. For the Earth- 
Sun system, at mid-transit, refraction prevents the remote observer's access to the lower 
~12-14 km of the atmosphere and, thus, also to the bulk of the spectroscopically-active 
atmospheric gases. We demonstrate that the effective optical radius of the Earth in 
transit is modulated by refraction and varies by ~12 km from mid-transit to 2nd con- 
tact. The refractive nature of atmospheres, a property which is rarely accounted for in 
published investigations, will pose additional challenges to the characterization of Earth- 
like extrasolar planets. Refraction may have a lesser impact for Earth- like extrasolar 
planets within the habitable zone of some M-type stars. 
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Introduction 



The contemplation of the Moon during a lunar eclipse reveals the dimming, and subsequent 
brightening, of the lunar disk as the satellite enters, and then exits, the shadow cast by the Sun- 
illuminated Earth. Ground-based observe rs have for decades constructed lightcurve s of the sunlight 
reflected from the eclipsed Moon (e.g., iLinkI Il962l : iGarcia Mufioz &: Pallg l201ll . and references 
therein). The experiment is relevant to the investigation of the terrestrial atmosphere because, 
for a known position of the Moon within the Earth's shadow, the lightcurves are int erpretable 
in terms of the optical properties of the atmosphere (lUgolnikov &: Maslovi 120061 . 120081 ). Beyond 
traditional photometric measurements, two spectroscopic observations of the August 2008 lunar 
eclipse have demonstrated the richness of the spectrum of sunlight reflected from t he eclipsed jMoon 
and shown that it cqntaiii s the signature of the principal atmospheric constituents (|Palle et al.ll2009l : 
Vidal-Madiar et al.l[201ol ^. 



Within the sample of known extrasolar planets, the transiting ones form a subgroup particularly 
apt for characterization. During transit, the planet partially blocks the disk of its host star, causing 
a drop in the amount of stellar light that arrives to the observer. The apparent stellar dimming is 
related to the sizes of the core and -when there is one- the gaseous envelope of the planet. The 
technique of in-transit transmission spectroscopy, which relies on comparing the apparent stellar 
dimming at selected wavelengths, has led to the detection of atoms and molecules such as Na, H, 
C, O, CO, CO2, HoO an d CH4 in the atmospheres of a few giant close- i n ext r asolar planets (e.g . 



Charbonneau et al 



Sing et al 



2009 



2002 



Vidal-Madiar et al.l l2003l . 12004 : iTinetti et al.l l2007l : iDesert et al.l l2009l : 

I 



Swain et al.ll2009l ). The same technique, app lied to GJ 1214b, is helping elucidate 

might contain abundant water in its 



whether this so-called super-E arth (iCharbonneau et al 
atmosphere ()Berta et al.ll2012l ). 



200£ 



Lunar eclipses and exoplanetary transits are related phenomena. In both instances a fraction of 
the light collected by the observer is stellar light that has passed through a section of the planet's 
limb and, thus, carries the signature of the planet's gaseous envelope. The photon trajectories, 
whether direct or deflected in scattering collisions with the atmospheric constituents, determine 
the signature imprinted on the collected stellar light. Throughout the text, we use the term (lunar) 
eclipse to refer to the alignment of the Sun, Earth and Moon, and transit to refer to the Earth's 
passage in front of the Sun as observed from a remote vantage point. 

The classical theory of lunar eclipses assumes that refr acted sunli ght , rather than the scattered 
compo nent, determines the brightness of the eclipsed Moon ( LinklE962 ). However, Garcia Mufioz et al 
(j201ll ) have shown that both components may become comparabl y intense under conditions of el 



evated aerosol loadiiig. W i th a few exceptions or brief references (ISeager &: Sasselovl I2OOOI: 



2001 



Hubbard et al. 



Sidis h Sari 



2001 



Hui fc SeageJ2002l : lEhrenreich et allboOfil : IXaltenegger fc TraublboOfll : 



Brown 



2OIOI ). the effect of refraction on the lightcurves of transiting extrasolar planets has 



largely been ignored. The omission is likely due to the minor contribution of refraction to the 
lightcurves of the close-in giant planets that constitute most of the extrasolar planets discovered to 
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date. However, and as shown below, refraction has a non-neghgible impact for smaher planets far 
out from their host stars. 



As the newly discovered extrasolar planets steadily approach terrestrial sizes (e.g . ICharbonneau et al 



2009; 


Leger et al. 


2009; 


Mavor et al. 


2009: 



tion is being devoted to the prediction of the in-transit signature of Earth-like extrasolar planets. 
In preparation for future observations, a number of recent works have investigate d the potential of 



in-transit transmission spectroscopy for the characterization of such targets (e.g. lEhrenreich et al 



20061 : iKaltenegger k TraublboogI : IPalle et ahlboill : iRauer et ablboilh . In that context, it has been 
suggested that lunar eclipse spectra may help t est the strategies that will eventually lead to the 
characterization of Earth- like extrasolar planets (jPalle et al.ll2009l ; IVidal-Madiar et al.ll20ld ). 



Our paper investigates the transport of stellar light through the Earth's atmosphere in lunar 
eclipses and exoplanetary transits. Our main goal is to clarify the similarities and dissimilarities 
in the signal collected by terrestrial observers in both types of event, emphasizing the impact of 
refraction. The paper is structured in two main blocks. In the first block, ^}2j-^|3l we introduce the 
lunar eclipse theory and investigate a set of previously unpublished near-infrared spectra of the 
August 2008 lunar eclipse. In the second block, ^ we take the lunar eclipse theory to the limit 
of a remotely distant observer and address the comparison between lunar eclipses and transits. In 
we also discuss the conclusions drawn by some prior works in relation with the theory of 
lunar eclipses and exoplanetary transits. Finally, ^summarizes the main conclusions. To the best 
of our knowledge, this is the first work that investigates the effect of refraction on the in-transit 
spectral signature of Earth-like extrasolar planets. 



The data. The lunar eclipse theory 



The spectra investigated in ^were o btained on 16 August 2 008 with the LIRIS spectrograph 
at the William Herschel Telescope, WHT (jManchado et al.lll998l ). located at the Observatorio del 
Roque de los Muchachos (ORM), La Palma, Spain. The observations were made alternating over 
two wavelength ranges with two different grisms: zj, from 0.9 to 1.5 /im, and /i/c, from 1.4 to 2.4 
//m. A total of 11 spectra of the Moon in umbra were collected, 6 in zj and 5 in hk. The umbra 
observations ran from 20:54 until 22:16UT and probed different phases of the Ea rth's inner shadow 
projected on the lunar disk. The data reduction and telluric correction follow iPalle et al.l (120091 ) 
and will not be discussed further. The FITS raw datafiles are publicly avail able from the ING 
Archive^ For their published spectrum from 1 to 2.4 //m, Palle et al. ( 2009 ) merged the three 
spectra closest to greatest eclipse in each wavelength range. We will now focus on the individual 
umbra spectra in order to investigate their evolution over time and distance from greatest eclipse. 

For the interpretation of the observations, we produced model simulations of the eclipse spec- 



^http://casu. ast.cam.ac.uk/casuadc/archives/ingarch 
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trum. The simulations involved solving the radiative transport problem of sunlight in the Earth's 



atmosphere according to the formulation laid out by iGarcia Muhoz &: Pallg (|201ll ) . The formu- 
lation assumes that the observer is located on the Moon's surface, which allows us to ignore the 
reflection of sunlight at the Moon and the atmospheric extinction above the Earth-bound observer. 

Generally, we integrate the direct and scattered components of radiance at the lunar observer's 
site over a solid angle di} that includes both the solar and planetary disks. Denoting by L(xo, sq) 
the radiance at the observer's site xq in the incident direction sq, the irradiance or flux of sunlight 
into an elementary surface oriented according to is obtained from: 

F(xo) = / L(xo,so)so • nxodfi(so). (1) 

Jan 

For the direct component, Beer-Lambert's law is integrated on the refraction-bent trajectories 
of all possible lines of sight connecting the lunar observer's site with the solar disk. At the solar 
disk, the Sun's emission radiance is scaled according to a specified limb-darkening function. For 
the scattered component, the solar photon trajectories are simulated by a Monte Carlo algorithm. 
The formulation is general enough to investigate both the umbra and penumbra phases of a lunar 
eclipse and, as discussed below, to predict the in-transit signature of extrasolar planets. 

Figure ([1]) sketches how the refraction-bent rays are traced sunwards from the observer's site 
and, in turn, how Eq. ([1]) is evaluated for the direct sunlight component. In our numerical imple- 
mentation of the integral, one ray is traced for each specified discrete element of sq • nxoC?^(so)- 
(We throughout refer to rays, although it would be more appropriate to refer to ray bundles.) 
Since the integration is carri ed out at x p, we do not have to correct for the optical phenomenon 



of attenuation by refraction (jLinklll969l ). Attenuation by refraction is simply the reduction in the 
cross section (and therefore in the associated irradiance) of a ray bundle departing from the solar 
disk and reaching the observer's site after crossing the Earth's atmosphere. In other formulations 
that use the solar plane for integrating Eq. ([1]), this phenomenon must be explicitly introduced as 
a factor wei ghing the size of the emitting solar disk parcel and its refracted image at the planet's 



terminator (lLinklll969l ). This effect would be especially important for rays crossing the lower layers 



of the atmosphere. Those are also the rays more strongly affected by gas and aerosol extinction. 

The formulation holds valid for any combination of e and do , where e is the geocentric angular 
distance from the area of lunar disk being probed to the geometrical umbra centre and do is the 
distance from the observer's site on the Moon to the Earth's centre. In the do—^oo limit, the 
formulation becomes relevant to the investigation of transiting extrasolar planets and e becomes 
a measure of the orbital phase. Numerically, we set do—^oo by assuming a sufficiently large do- 
Special care is taken to ensure that all the formula occurring in the numerical implementation adopt 
the correct forms in that limit case. 
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The eclipse in progress 



We start by investigating the near-infrared spectra obtained at the WHT, looking for changes 
in their structure during the umbra phase of the echpse. The exercise is useful to further validate 
our model of lunar eclipses and to show how the depth of molecular bands in the spectra varies as 
the eclipse progresses. 

The nominal atmosphere in the simulations takes the tem perature, O3 and Hg O (vapor) pro- 



files from the FSCATM subarctic summer model atmosphere (jGallerv et alJIlQSsI). and constant 
volum e mixing ratios of 0.2094, S.SSSxlO""^ and 1.779x10"^ for O2, CO2 and CH4 (jHansen Sato 



2004 ) , respectively. The vertical profiles of temperature and the volume mixing ratios are graphed 
in Fig. ([2]) . The density at all altitudes was determined by integration of the hydrostatic balance 
equation. For the Collision Induced Absorption (CIA) bands of oxygen at 1.06 and 1.27 ^um, we 
proceeded as follows. For the 02-02 CIA band at 1. 06 /xm, we adopted t he bi nary cross sections 
measured at 230 K in a 75/25% mixture of O2/N2 bv lSmith k NewnhamI (|20od ). For t he O^-O^ + 



O9 N0 CIA band at 1.27 //m, we used the binary cross sections measured at 253 K by lMate et al 



( I999I ). Both sets were rescaled to have peak values in air of 1.7x10-^^^x0.21=3.57x10-^^ and 
1.7x10-^^x0.21x4=1.43x10-^^ cm^ molec-^ at 1.06 and 1.27 //m, respectively. The selected 
band shapes and scaling factors are consistent with the conclusions drawn in a recent investiga- 
tion of solar occultation data obt ained with the SCIAMACHY spectrometer aboard ENVISAT 
(|Garcfa Muhoz &: Bramstedtll2012l ). We evaluated the optical opacity for both CIA bands through 
dr = (Tcia[02] [X](is, where ds is the differential integration path, [•] stands for number density, X 
is the O2 collision partner in the CIA band (O2 at 1.06 //m and O2 + N2 at 1.27 //m) and aciA are 
the binary cross sections. We set an opaque cloud layer with cloud to ps at 6 km and an aerosol ex- 
tinction profile at 1.02 /im about 4 times the September 2008 profile of lSioris et al.l ()2010l ). Further, 
we assumed that the aerosol optical properties are wavelength-independent ov er the explored range 



of_wavelengths. The nominal atmosphere is consistent with the findings by iGarcia Muhoz et al 



(2011 



from optical data of the same eclipse that point to a heavy aerosol loading following the 
eruption of the Kasatochi volcano one week ahead of the eclipse. Scattered sunlight is negligible 
longwards of 1 fim and therefore omitted. The simulated spectra were convolved with a Gaussian 
line shape at a resolving power of about 1,000. 

A critical aspect in the comparison between observations and simulations is the solar elevation 
angle, e, or distance from the geometrical umbra centre to the projection of the instrument's 
collecting element on the lunar disk. No images were taken of the slit projected on the Moon and 
e had to be estimated indirectly. For that purpose, we followed two different methods. Method 1) 
estimates e from the telescope pointing coordinates recorded in the FITS raw datafiles. The Sun, 
Earth and Moon coordinates, the coordinates a.t the ORM and the angular size of the lunar disk 
were taken from the JPL HORIZONS service (jGiorgini et al.lll996l ). In method 2), we produced 
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geocentric images of the lunar disk with the GeoViz softwarqj. From the Moon's centre, we displaced 
the slit south in equatorial declination on the image until only half of the 4.2-arcmin slit rested on 
the lunar disk. From the image, we recorded the coordinates of the mid-point projection of the slit 
on the lunar disk. Then, e was calculated as the geocentric angle from that point to the umbra 
centre. 

Figure ([3]) shows the set of 11 umbra spectra together with our model simulations. We arranged 
the zj and hk spectra in pairs according to the times of observation. Each pair contains data 
collected within an interval of no more than 6 minutes. To first approximation, each pair can be 
seen as an uninterrupted spectrum from 0.9 to 2.4 /im. For the model simulations, we adjusted 
the e parameter to ensure an optimal match of the eclipse data. In some cases we modified the 
H2O content by an amount that was always less than 10% of the nominal profile. The text in 
the graphs gives the time at mid-exposure and airmass for the Moon-to-telescope optical path, 
along with the estimated (ei and 62) and adjusted (eadj) solar elevation angles. The quality of 
each match was judged by giving special weight to the O2 and CIA bands in zj, and the CO2 and 
CII4 bands in hk. Absorption in both CIA bands remains linear throughout the eclipse. Indeed, 
we note the suitability of the CIA band at 1.06 ^um, which occurs over a region only moderately 
affected by other molecular bands, for monitoring the Moon's progress across the umbra. Overall, 
the agreement between the eclipse data and the simulations is consistently good over the entire 
spectral range. Exceptions are the regions of strong water absorption at 1.35-1.45 and 1.8-2.1 /im, 
where the telluric correction introduces obvious artifacts. A more flexible implementation of aerosol 
extinction in the simulations would have surely improved the match between the observational and 
synthetic spectra, but it was judged that such an effort was not critical for the present purpose. 

The minimum Cadj occurs at 21:12-21:16 UT, in fair agreement with the prediction of geomet- 
rical greatest eclipse for 21:10:06Ut|£] The comparison between the estimated and adjusted angles 
indicates that ei and 62 are systematically larger than Cadj- The differences are typically ~0.05° 
except in the 20:54-20:59 UT spectrum, in which case it reaches ~0.1°. This spectrum was acquired 
at a rather high airmass, and refraction in the Moon-to-telescope optical path may thus play a role. 
The 4.2-arcmin slit spans ~l/15-th of the lunar disk diameter. The sunlight that enters the slit 
arrives from a range of distances to the umbra centre that may lead to differences in e of ~0.04°. 
Further, the non-uniform albedo of the Moon along the slit may move the effective slit centre away 
from its geometrical location. It is thus difficult to define an equivalent solar elevation angle from 
purely geometrical considerations. 

In addition to the above arguments, one may also expect that local features in the atmosphere 
intercepted by the refracted sunlight rays will have an impact on Cadj- To explore this possibility, 
we conducted a few simulations in modified conditions of clouds and aerosols. As a rule, setting 



^New Horizons Geoviz by H. Throop, http: / /soc.boulder. swri.edu / nhgy/] 
^Eclipse Predictions by F. Espenak, |http: / /eclipse. gsfc.nasa.gov/eclipse. html 
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the cloud tops below the nominal 6-km level leads to deeper absorption bands than in nominal 
conditions. If, for instance, the cloud tops are set at 3 km, the 21:12-21:16UT eclipse spectra are 
optimally reproduced with e'^0.29°, rather than 0.27°. The effect of low clouds on eadj diminishes 
progressively as e increases. The nominal aerosol profile extends well up to ~ 16-1 7 km and results in 
sunlight extinction over a broa d range of alt i tudes . We verified that reducing the aerosol extinction 



to the September 2008 levels of lSioris et al.l (j20ld ) has a minor impact on the spectrum's structure 
near greatest eclipse. 

Figure dH top) shows the 21:12-21:16UT eclipse spectrum and the corresponding simulation. 
In Fig. (jH middle), the color lines are the contributions from H2O, O2, CO2, CII4 and the CIA 
bands of oxygen to the simulation. For e~0.3°, the sunlight rays that reach the area of lunar disk 
bein g probed have their closest approach to the Earth's surface at altitudes in the range ~2-ll 



km (jGarcia Mufioz fc Pallell201ll ). The range is effectively narrower because the lower altitudes are 
optically thick and, in addition, we assumed opaque clouds below 6 km. The spectrum of Fig. ([H 
top) is indeed well approximated by a limb- viewing transmission spectrum of the atmosphere as seen 
from a tangent altitude of ~10 km. For comparison, Fig. dH bottom) shows two solar occultation 
spectra measured with the SCIAMACHY spectrometer for tangent altitudes of ~7.1 and 9.5 km. 
The resemblance of the lunar eclipse spectrum to the SCIAMACHY spectra is apparent. 

For e~0.5° (still in the umbra), the range of altitudes of closest approach is ~4-16 km, and 
for e~0.7° (at the umbra/penumbra edge) it is ~8~65 km. The shift and widening of the range of 
altitudes for closest approach has notable consequences on the spectrum's structure. For molecules 
whose densities decay monotonically with altitude, increasing e means the dilution in the spectrum 
of the molecular signature. This is due, first, to a diminished contribution from the atmosphere's 
optically thicker layers and, second, to a larger amount of light passing unattenuated through the 
upper altitudes. This effect is especially apparent for H2O and the 02-02 + 02-N2 collision complex 
because their densities drop at a faster rate than those of (nearly) well-mixed molecules such as 
O2, CO2 and CH4. The simulations of the monomer and CIA bands of oxygen at 1.27 /im in Fig. 
([5]) illustrate this trend. For e=0.3°, the 02:CIA equivalent widths are in a ratio of 1:3.3, whereas 
for e=0.7°, the ratio is only 1:0.9. Thus, as the observer probes regions of the lunar disk closer 
to the penumbra, the relative contributions of II2O and the collision complex with respect to the 
other molecules tend to decrease. 



3.1. The impact of the upper atmosphere on umbra spectra 



The umbra spectrum published bv lPalle et al.l (1200911 with data of the 16 August 2008 event is 
representative of conditions near greatest eclipse (jGarcia Muhoz et al.ll201ll ) and the atmosphere at 
mid-to- north latitudes over the Atlantic. At the relevant e~0.3°, the absorption features are formed 
in refracted sunlight ray trajectories that approach the Earth's surface at minimum distances <11 
km. The occurrence of strong H2O and CIA band features confirms the importance of the lower 
altitudes. 
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Palle et alj (|2009l ) mention the identification of the Na I neutral atom at optical wavelengths 
in their published spectrum. The eclipse was contemporary with the meteor shower of the Per- 
seids, which means that abnorma lly elevated amounts of that metal in the m esosphere and lower 
thermosphere might be expected (|Planell2003l ). Following iFussen et alJ (|2004l ). we estimate that in 
globally-averaged conditions the optical thickness at the resolution of the umbra spectrum for the 
Na I doublet at 0.589 fim from a tangent altitude of 11 km is ~2xl0^'^. The sodium peak density 
is known to increase sporadically by factors of up to 10 in layers a few kilometer thick (which 
entails an increase in the integrated column by a factor of a fe w), sometimes over hor i zontal spans 
of hundreds of kilometers, especially during meteor showers (jMoussaoui et al.l l20ld : iPland l2003l : 



Dou et al.ll2009l ). Even then, such a weak signature makes the detection of the Na I neutral atom 



challenging for a spectrum of moderate spectral res o lution . Thus, our analysis indicates that the 
identification of the Na I doublet in the I Palle et al.l (|2009l ) spectrum is questionable and requires 
additional confirmation. 



4. The Earth in transit 

4.1. The Earth-Sun system at mid-transit 

Everyday, sunsets offer real-life demonstrations of the refraction of sunlight in the atmosphere. 
Looking from the ground over the horizon, the setting Sun appears about one solar diameter 
above its true elevation. Sunlight rays crossing the atmosphere and having their closest ap- 
proach to the Earth's surface at a tangent altitude /tf ^ p are approximately refr acted by an angle 



arefr~(ra(/itan)-l)x(27ri2p/i/)i/2 (|Baum k Codelll953l : ISeager &: Sasselovll2000l ). For a terrestrial 
radius i?p~6377 km, a constant scale height ff~8 km and a typical refractivity at sea level in the 
optical n(/itan=0) — l~2.8x 10~^, the end-to-end refracted angle estimated from the above equation 
is Orefr'^l-l^". Thus, over a distance of 1 AU, a ray incident on the planet that grazes the Earth's 
surface is deflected by ~4.26 solar radii. Tracing the ray trajectories from a remote observer's site 
back towards the solar disk, one concludes that the mid-transit solar image at the observer's site is 
formed by sunlight rays having their closest approach to the planet's surface at densities ~l/4.26 
or less than at sea level. Further, our estimate suggests that, at mid-transit, refraction prevents the 
remote observer's access to altitudes within a Zrefr-sized ring above the planet's surface. The exact 
size of this refraction-exclusion ring can only be determined by integrating the ray trajectories 
through the atmosphere after prescribing the altitude-dependent refractivity of the atmospheric 
gas. Assuming, for the purpose of obtaining a first estimate, that the atmosphere is isothermal and 
that the refractivity decays with altitude as exp (—z/H), one obtains that Zi.efr~-f^lii4.26'^11.6 km. 
This approximate result is reached by equating the deflected distance of rays tangent at an altitude 
ZreiT and the solar radius, i.e. Rq = (n(/itan = 0) — 1) exp (— Zi.efr/-ff)(27ri?p/i7)^/^a0. Here, Rq 
and a© stand for solar radius and the planet's orbital distance, respectively. 

Our model simulations with the formulation of ^in the do^oo limit provide refined magni- 
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tudes of the above estimates for Zrcfr- For the nominal atmosphere described above, a sunhght ray 
grazing the Earth's surface is refracted by ~1.08° and only altitudes above 2;refr'^13.2 km contribute 
to the mid-transit signal measured at the observer's site. Thus, the mid-transit spectrum measured 
by the remote observer will lack the signature of the atmospheric layers containing the bulk of the 
spectroscopically-active gases. The value of Zj-efr depends on the fluctuations in the atmosphere's 
density profile and water content through the refractivity of the gas. A sen sitivity analys i s per - 
turbing those two parameters within the values in the model atmospheres of Gallery et al. (jl983l ) 
shows that Zrefr typically lies between 12 and 14 km. 

Figure ^ displays in black the mid-transit spectrum of the Earth in two levels of approx- 
imation, namely with and without refraction, and three atmospheric scenarios. The scenarios 
represent: (1) a Rayleigh atmosphere, free of clouds and aerosols; (2) an atmosphere with cloud 
to ps at 2 km, the bac kground aerosol extinction profile for September 2007 at about 1 //m published 
bv lSioris et al.l (j20ld ) multiplied by (1.02/A [/im])^-^ at other wavelengths; and, (3) the aerosol-rich, 
cloudy atmosphere described in fj3l The aerosol loading in scenario (2) might be seen as a plausible 



representation of globally-averaged conditions (jHayashida &: Horikawall200ll ). The aerosol loading 



in scenario (3) is abnormally elevated for typi cal Earth conditions but ni ay be representative of the 



atmosphere after a major volcanic eruption (jGarcia Muhoz et al.ll201ll ). We note, however, that 



in our own Solar System, Mars undergoes episodic events of abnormally elevated aerosol amounts 
being transported across the globe and thus affecting the aerosol loading on global scales, and that 
Venus is enshrouded by a complex system of clouds and haze layers up to pressures well below 1 
atm. 



We use the equivalent height, heq, defined according to: 



1 



Rp + h, 



cq 



(2) 



to provide a measure of the atmosphere's thickness opaque to the incident sunlight. As a consistency 
check, the refractionless calculations were te sted against the more usual formulation for the in- 
transit stellar dimming ( Hubbard et al. 200ll ): 



Fp _ 2vr/^° exp(-r(rb))r5drb 



^0 



with: 



rin) 



R, 



d{r/n), 



(3) 



(4) 



Rp and 7(r) is the optical extinc- 



where r(rf,) is the optical thickness at an impact altitude 
tion coefficient at altitude r — Rp above the planet's surface. The agreement between the two 
implementations proved to be excellent. It must be mentioned that Eq. ^ was integrated in t, 
=cosh~^ (r/rh), rather than in r/r^, to avoid the singularity in the denominator of the integrand. 

Figure Q clearly shows how refraction removes much of the mid-transit spectrum's structure. 
An exception to that is the Chappuis band of ozone, which absorbs noticeably in all cases shortwards 
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of ~0.8 /im. Terrestrial ozone densities peak in the stratosphere, and ozone in the refraction- 
exclusion ring contributes in a minor way to the disk-integrated signature. The features longwards 
of 1 fj,m are readily identifiable by comparison with Fig. Shortwards of 1 fj,m, the sharp 

features at 0.76, 0.69 and 0.63 /im are the X{v"= 0)^b{v'= 0, 1, 2) bands of O2, respectively. The 
remaining structure in that region is largely due to H2O. With respect to the nearby continuum, the 
equivalent heights of the discrete features in the spectra of the refractive atmosphere are always less 
than ~10 km in the 0.4-2.5 /im spectral range. In going from atmospheric scenarios (1) to (3) it is 
seen that refraction becomes less important as aerosol extinction near the Zrpfr boundary b e come s 



do minant. The visual inspection of the transmission spectra published bv lEhrenreich et al.l (|2006l ) 



and lKaltenegger &: Traubl (j2009l ) for an Earth twin shows that their spectra are generally consistent 



with our refractionless calculations. 

Were the Sun and Earth observed at mid-transit from a remote distance, the transmitted 
sunlight reaching the remote observer's site would be determined by sunlight photons crossing 
the Earth's atmosphere at altitudes above ~12-14 km. As seen in Fig. ([6]), ignoring the refractive 
nature of the Earth's atmosphere would mean a significant overestimation of most spectral features. 
The overestimation is more obvious when the model atmosphere is assumed free of clouds and 
aerosols since airborne particles act as natural barriers to the transmission of light through the 
lower altitudes. Because meteorological activity is largely confined to the troposphere, the planet's 
meteorology will have limited impact on the in-transit transmission spectrum. 

The sunlight rays passing closer to the planet's surface at mid-transit originate from the near- 
limb region of the opposite solar hemisphere. In order to assess the impact of solar limb darkening, 
we repeated the calculation for the three scenarios described a bove including a linear limb-darkening 



law for the solar brightness, U = 1 — ui{l — Hq), and ui=0.6 (!Gimeneal2006l ). To proceed with the 
comparison, we redefined h^q in the way: 



Fq 3-ui \ Rq 

to account for the (3— ni)/3 factor in the limb-darkened solar irradiance, Fq. Our model simulations 
indicate that heq differs by less than ~1 km in the cases with ui= and 0.6. Thus, solar limb 
darkening affects in a negligible way the spectrum's structure at mid-transit. 

As in ^ we calculated the equivalent widths of the monomer and CIA bands of oxygen at 1.27 
/im. For the refractive calculations of scenarios (l)-(3), we obtained that the 02:CIA equivalent 
widths are in ratios of 1:0.70, 1:0.64 and 1:0.47, respectively. The ratios indicate that the distinct 
identification of the CIA band becomes more and more difficult as the atmospheric opacity of 
aerosols and clouds increases. 



Garcia Munoz &: Pallg ((20111) estimated that the sunlight scattered at the Earth's terminator 
towards the eclipsed Moon could amount to up to ~10~^-10~^ of the solar irradiance depending 
on the forward-scattering efficiency of the airborne aerosols. During the eclipse, the magnitude of 
scattered sunlight depends on the ratio of the solid angles subtended by the Earth, ~l/dQ, and the 
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Sun, ~l/((io + fl©)^) from the lunar observer's site on the Moon. In the lunar eclipse, we used for 
do and gq, 382665.9 km and 1 Astronomical Unit, respectively. Thus, the amount of stellar light 
scattered towards the remote observer during the transit is ~(10~''-10~^) x (do + a0)^/(^o~ 10~^^- 
10~^^, which is much less than the ~10~^ dimming caused by the planet's core or the 10~^-10~^ 
dimming attributable to the strongest atmospheric features. 

Next, we address under what conditions a lunar eclipse spectrum may be representative of the 
in-transit transmission spectrum of the remotely-observed Earth-Sun system. Since at mid-transit 
only altitudes above ~12-14 km are probed, the eclipse phases matching mor e closely that condition 



entai l solar elevation angles e>0.9°, which are well within the penumbra (jGarcia Muhoz &: Palle 



20 111 ). Thus, spectra from the penumbra, rather than from the umbra, will b e better representatives 



of the in-transit transmission spectrum for the Earth-Sun system. Figure 3 in lGarcia Mufioz &: Palle 



(|201ll ) shows an eclipse spectrum at different phases. The comparison of that figure and our Fig. 
([6]) confirms that the eclipse spectrum tends to resemble the mid-transit transmission spectrum as 
the eclipse approaches the umbra/penumbra edge. 



4.2. The Earth-Sun system near ingress/egress 



We have so far focused on the spectrum at mid-transit, noting the effect of refraction on the 
equivalent height of the Earth's atmosphere with respect to the situation for a virtual non-refractive 
atmosphere. The refraction-bending of sunlight rays does, however, introduce phase-dependent 
effects on the intensity and spectral structure of the transmitted sunlight that are worth exploring. 

Out of transit, refraction causes the brightening of the planet's hemisphere more distant from 
the solar disk. This effect has been reported for Venus transits and causes a transient halo at 
the planet's limb that is readily discernible from Earth with the aid of a small telescope between 
internal and external contacts. It is attributed to Lomonosov back in the 18th century the correct 
interpr etation of t his halo, a fact that, in turn, meant the first evidence for an atmosphere at the 
planet (|Linklll959l ). The last Venus transit occurred in June 2004 and provided numerous images 
of the halo t hat confirmed the pole-to-equator asymni etry in the optical properties of the planet's 
atmosphere (jPasachoff et al.l l201ll : iTanga et al.l l2012l ) . The contribution of the halo to the net 
brightness of the planet-star system is tiny for Venus but might amount to a detec table magnitude 
for transiting Jovian extrasolar planets on long-duration orbits (jSidis &: Sarill2010l ). 



The halo is a visual representation of the enhancement in the amount of refracted sunlight at the 
planet's outer hemisphere. This enhancement occurs together with a diminishment on the opposite 
hemisphere. Figure ^ sketches the situation. Sunlight rays passing through the outer hemisphere 
may be refracted by an angle of up to ^2Rq/ a^, whereas on the inner hemisphere refraction angles 
of more than ~2i?p/a0 do not contribute to the signal collected by the distant observer. Different 
refraction angles also mean different depths of the sunlight rays into the atmosphere and, in turn, 
a refraction ring whose size is larger on the inner hemisphere than on the outer one. 
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The upper curves in Fig. ^ represent the equivalent heights of the atmosphere from mid- 
transit (e=0°) to internal contact (e~0.2641°). We omitted limb darkening to emphasize the purely 
refractive effects, an assumption that enables us to use Eq. ([2]) as a valid definition for the effective 
size of the atmosphere. The calculations assumed the Rayleigh atmosphere described above. The 
main conclusion to draw from that set of curves is that the Earth's optical size varies with the phase 
of the transit as a result of the refraction of sunlight at the planet's limb. The Earth appears larger 
near internal contact than at mid-transit because in the proximity of the solar limb the optical 
enlargement of the inner hemisphere dominates over the shrinking of the outer one. The overall 
change in Fp/FQ from mid-transit to internal contact amounts to about 3xl0~^ for a variation 
in the equivalent height of 12 km. The detection of such a dimming in the stellar light appears 
extremely challenging, especially after considering that most of the drop occurs near ingress/egress 
where limb darkening effects on the planet-star lightcurve are more prominent. The upper set of 
curves in Fig. ([8|) also reveals that all molecular band features tend to weaken for phases that move 
away from mid-transit. 

The lower set of curves in Fig. ([8]) represent a measure of the halo brightness from external 
contact (ew0.2690°) to phases well out of the eclipse. We use an equivalent halo thickness defined 
through: 

-Fp _ , 27ri?p/lhalo /p-N 

that gives the thickness of a ring around the planet of the same brightness as the Sun for /iQ=l. The 
effective halo thickness is about 10 km at external contact and decays as the planet moves farther 
out of transit. The net brightening at external contact is about 3x10^"^ and, again, probably too 
faint for detection with the currently existing technology. 



4.3. An Earth in transit in a different planetary system 

We have seen that at mid-transit refraction imposes a ring around the Earth's terminator that 
deflects some line of sights from courses reaching the solar disk or, equivalently, some sunlight 
rays incident on the planet from courses reaching the remote observer. The size of this refraction- 
exclusion ring depends on the refractive properties of the atmosphere but also on the stellar radius 
and the planet's orbital distance. In standard conditions of temperature and pressure, the N2 
and O2 refractivities are about 1.6 times less than the CO2 refractivity but about twice the H2 
refractivity. For a virtual Earth twin whose atmospheric properties remained unchanged when 
placed in a planetary system other than ours, the ring will shrink if the planet is closer than 1 
AU to the star {as < Uq) or if the host star is larger than one solar radius {Rg > Rq)- The 
inequality arefras/^s> 1 for the ratio of the refraction-deflected distance, arefr^s (with ai.efr~l°)) 
and the stellar radius establishes the approximate geometrical condition for which a refraction- 
exclusion ring exists at sub-atmospheric pressures. Speculations on possible values for arefvCi's/Rs 
find their justification in the wide variety of known planetary systems and the even wider variety 
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that is likely to be discovered. As a matter of fact, Fig. ^ applies to a broader range of ag/aQ 
and Rs/Rq parameters than first stated. The lower synthetic curve represents approximately 
the condition as/aQxRQ/ Rs<l and small enough so that refraction has a minor impact on the 
in-transit transmission spectrum. In turn, the upper synthetic curve represents the condition 
ag/aQ^RQ/ Rg'^l. Those two curves bracket all possible conditions for the intermediate values of 
as/aQy.RQ/Rs. 

Oxygen is a potential biomarker and therefore constitutes a particularly appealling target for 
future searches. From our Fig. ([6]), it is seen that the O2 X(0)— )-6(0) band at 0.76 /im may be 
the best candidate for an oxygen search, although the CIA bands at 1.06 and 1.27 ;um (the latter 
blended with the simultaneous monomer band) are also good candidates when the lowermost layers 
of the in-transit extrasolar planet are probed. 

Although refraction will generally pose an additional challenge to the characterization of Earth- 
sized planets, the added difficulty may be minor in some configurations of interest. That is indeed 
the case for an Earth twin orbiting within the habitable zone (HZ) of some M-type stars, a planet- 
star combination that yields relatively favourable areal ratios. The HZ is the distance around 
a star in which a planet's water might be expected to exist i n liqu id state. Taking for the HZ 



the reference orbital distance quoted by iKaltenegger &: Traubl ([20091), ag/aQ = {Ts/Tq)'^{Rs/Rq), 



where, following our earlier convention, subscripts s and denote stellar and solar magnitudes, 
respectively, one arrives at the condition, independent of Rs, (T^/Tq)^ < {Rq / uq) / Uj^cir for the 
stellar effective temperature below which refraction effects are minor. If we adopt arefr= 1.08° at 
the 1-atm pressure level, the condition becomes Ts/Tq<0.5. Thus, an Earth twin that orbited 
within the HZ of an M-type star of effective temperature about 2900 K would produce mid-transit 
spectra with a closer resemblance to the refractionless simulations in Fig. ([6|) than to the refractive 
simulation of the Earth-Sun system. Obviously, other intermediate situations would be possible for 
Ts values in between 2900 and 5800 K. 



Discu ssing detectability issu e s for s pecific telescope sizes in the way it is d one bvlEhrenreich et al 

2006 ) or IKaltenegger &: Traubl ([20091) is beyond the scope of our work. IKaltenegger fc Traubl 



2OO9I ) do recognize the potential impact of refraction on the lightcurves of Earth-sized planets. 



They do not, however, give the size of the refraction ring in their investigation of the in-transit 
transmission spectrum of an Earth twin orbiting M-type stars nor explain the way refraction is 
implemented. The authors mention an altitude of 6 km below which the terrestrial atmosphere is 
opaque to grazing sunlight rays but do not state how much of that is contributed by refraction. 
As pointed out above, the size of the refraction-exclusion ring depends on the refractive proper- 
ties of the atmos p here but also on both the stellar radius and the planet's orbital distance. Since 
Ehrenreich et al.l (120061 ) account in some cases for the presence of clouds at altitudes of up to 10 
km, the conclusions drawn from their cloudy models of an Earth twin might be appropriate for 
specific conditions in which refraction plays a significant role. 



Palle et al.l (j201ll ) discussed the in-transit detectability of the strongest absorption features 
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from 0.4 to 2.4 /im using scaled versions of the iPalle et al.l (|2009l ) umbra spectrum. Setting the 
noise terms to zero in Eq. 4 of iPahe et al.l (|201ll ). the equivalent height proposed in that work 
can be written as /icq = /itoa x (1 — TTp), where Htoa is a specified top-of-the-atmosphere height 
and Tp is the umbra spectrum normalized to one near 2.2 ^m. Our Fig. ( [6]) shows in r ed th e 
resulting equivalent heights for the specific case /itoa =40 km investigated in iPalle et al.l (j201ll ). 
The comparison with the upper (black) synthetic curve in the top panel demonstrates that for an 
Earth-like extrasolar planet lacking aerosols and cloud s and close enough to a solar-type star so that 
refraction effects are minor, the approach followed bv lPalle et al.l (|201lh overestimates moderately 
th e model prediction s for a Rayleigh atmosphere. This is because the umbra spectrum published 
bv lPalle et al.l (|2009l ) results from probing preferentially the lowermost altitudes of the atmosphere 
during the eclipse. In all other cases displayed in Fig. ([6]), a nd more especiall y when refraction 
effects play a role, the equivalent heights for /itoa =40 km in lPalle et al.l (|201ll ) are larger by up 
to an order of magnitude than the current model predictions. It is apparent that fainter features 
will impose tougher requirements on their detectability. 



5. Concluding remarks 

We investigated lunar eclipses of the Earth and the transits of Earth-like extrasolar planets in a 
comparative manner. The Moon in umbra receives unscattered sunlight that has been transmitted 
through a section of the terrestrial atmosphere dictated by refraction. The diffuse sunlight that 
arrives at the eclipsed Moon is scattered from the entire terminator. During the transit of the 
Earth-Sun system, at mid-event, the refracted sunlight is transmitted through an annular ring that 
fully encloses the planet. The refraction-bending of sunlight rays prevents sunlight rays passing 
through a refraction-exclusion ring above the planet's surface from reaching the remote observer. 
For the Earth-Sun system, at mid-transit, this ring extends up to 12-14 km, thus blocking the 
observer's access to the denser atmospheric layers. It may be stated that umbra eclipses maximize 
the contribution of the lowermost atmosphere and, to some extent, of scattered sunlight, whereas 
during the transit of an Earth-like extrasolar planet the transmitted signal is dominated by altitudes 
above a refraction-exclusion ring and scattered sunlight becomes negligible. 

In the course of this work it has become apparent that a genuine transmission spectrum of the 
Earth in transit as observed from a remote distance is more difficult to obtain than what had been 
anticipated. Lunar eclipse spectra obtained during the penumbra are more representative of the 
in-transit spectrum than those obtained during the umbra. In the context of Earth-like extrasolar 
planets, the lunar eclipse spectra are helpful to fine tune synthetic spectra incorporating the fairly 
good knowledge of the atmosphere obtained over decades of Earth's atmospheric research. 

Refraction will affect the lightcurves of transiting Earth-like extrasolar planets. Its significance 
will depend on the refractive properties of the planet's atmosphere, the stellar size and the planet's 
orbital distance. All these factors must be considered in the future to predict and/or interpret the 
lightcurves of such planets. The role of refraction will not be critical for Earth twins orbiting within 
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the HZ of M-type stars with effective temperatures of about 2900 K and less. Finally, it is worth 
noting that recent lunar eclipse observations have contributed critically to bring up the importance 
of refraction in the future characterization of Earth-like extrasolar planets. 



The near-infrared spectra discussed in ^were obtained with the William Herschel Telescope 
at the Spanish Observatorio del Roque de los Muchachos of the Instituto de Astrofisica de Canarias. 
MRZO and PMR acknowledge support from projects AYA2010-21308-C03-02 and CGL2009-10641, 
respectively. Both MRZO and ELM acknowledge support from project AYA2011-30147-C03-03. 
ELM is grateful for financial support from RoPACS, a Marie Curie Initial Training Network funded 
by the European Commission's Seventh Framework Program, and b y the CONSOLIDE R-INGENIO 
GTC project. AGM acknowledges G. T. Eraser for providing the I Mate et al.l (jl999l ) binary cross 
sections and K. Bramstedt and the European Space Agency for the SCIAMACHY solar occultation 
spectra. 
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Fig. 1. — Sketch of the Earth-Sun system and the ray tracing scheme from the observer's site 
towards the Sun. Each ray is characterized by an incident direction so at the observer's site xq. 
The elementary surface for evaluating the irradiance at xq is oriented with its normal vector 
following the line that joins xq and the planet's centre. 
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Fig. 2. — Vertical profiles of temperature (dashed curve) and volume mixing ratios (solid and dotted 
curves) in the nominal atmosphere 
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Fig. 3. — In red, the eclipse spectra in zj (left) and hk (right). In black, the model simulations. 
The text inset shows the estimated e angles and the e angle adjusted in the model simulation to 
match optimally the observed spectra. AM is the airmass for the Moon-to-telescope optical path. 
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Fig. 4. — Top: Lunar eclipse spectrum at 21:12-21:16UT and model simulation. Middle: O2, 
02-X, H2O, CO2 and CH4 contributions to the model simulation. Bottom: SCIAMACHY solar 
occultation spectra measured from tangent altitudes of 7.1 and 9.5 km. 




1 .1 5 1 .20 1 .25 1 .30 1 .35 
Wavelength [/xm] 



1 .20 1 .25 1 .30 1 .35 1 .40 
Wavelength [/xm] 



Fig. 5. — The O2 (solid) and CIA (dotted) absorption bands at 1.27 /im for e angles of 0.3 and 0.7°. 
Far from the umbra centre, the CIA band loses relative importance with respect to the monomer 
band. 
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Fig. 6. — Equivalent heights for a refractive and a non-refractive atmosphere. Top panel: Rayleigh 
atmosphere. Middle panel: Globally-averaged conditions of clouds and aerosols. Bottom panel: 
Aerosol-rich, cloudy atmosphere as described in f|3j The red curve on each panel is the equivalent 



height as defined through the noise- free version of Eq. 4 in lPalle et al. (l201lh: ho„ = ^t qa x (1 — T 



with /iTOA=40 km and where Tp is the umbra spectrum published bv lPalle et al.l (j2009l ) normalized 
at its maximum near 2.2 fim. 
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Fig. 7. — The Earth-Sun system at mid-transit and at 2nd contact. In each case, the sunhght that 
traverses the Earth's Hmb originates from a different zone of the solar disk. At 2nd contact, the 
outer hemisphere sees the entire solar disk, whereas the inner hemisphere sees only a narrow zone 
close to the Sun's edge. 
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Fig. 8. — Equivalent heights for the Earth-Sun system in transit (upper set of curves) and halo 
heights for the system off transit (lower set of curves) at various phases (e). Limb darkening has 
been omitted to highlight the impact of refraction. For the in-transit curves, hgq is defined through 
Eq. ([2]), whereas for the off-transit curves, heq is defined through Eq. ([5]). The calculations in this 
figure are specific to a Rayleigh atmosphere. 



